[z,xc,yc]=exp_surface(-0.02+0j,[-0.04-0.1*pi*j 0+0.1*pi*j],100,[0:100]');
figure(2) ; plot(yc,z(100,:))
figure(1) ; mesh(yc,xc,z) 